Critical behavior of colloid-polymer mixtures in random porous media 
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We show that the critical behavior of a colloid-polymer mixture inside a random porous matrix 
of quenched hard spheres belongs to the universality class of the random-field Ising model. We 
also demonstrate that random-field effects in colloid-polymer mixtures are surprisingly strong. This 
makes these systems attractive candidates to study random-field behavior experimentally. 
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One longstanding problem in the phase behavior of flu- 
ids concerns the universality class of the liquid-gas tran- 
sition in random porous media. Conceivable universality 
classes are those of the pure Ising model, the (bond or 
site) diluted Ising model 0, or the random-field Ising 
model 1^ 9 ' It was suggested by de Gennes ^ that 
the universality class is that of the random-field Ising 
model. However, resolving the universality class exper- 
imentally is difficult, partly because well-characterized 
porous media are scarce. The prototype realization is 
silica aerogel, which has the disadvantage that the cou- 
pling between the porous medium and the fluid is weak. 
This is evident from the small shift in the critical tem- 
perature observed in these systems 0, 0| 

<5=|TM-Tp|/^p«6xlO-^ (1) 

where Tp is the critical temperature in the pure sys- 
tem, and Tm the critical temperature in the aerogel ma- 
trix. The universality class, consequently, could not be 
resolved in these experiments. 

In theoretical approaches, the porous medium is usu- 
ally modeled by an equilibrium configuration of fixed 
spheres. The fiuid particles are then allowed to migrate 
through the medium. Interestingly, these quenched- 
annealed systems display pronounced shifts in the critical 
temperature, typically 5>O.2@,[i0,EiEli3- Com- 
pared to aerogel, the coupling between the fluid and the 
porous medium is thus much stronger, making quenched- 
annealed fluids attractive model systems. 

Nevertheless, the universality class of quenched- 
annealed systems has not yet been determined. The 
most direct approach, which is to measure the critical 
exponents, is cumbersome because the critical exponents 
of the pure and diluted Ising model are rather similar, 
whereas for the random-field Ising model (RFIM), the 
critical exponents are not precisely known in any case. 
However, there is an additional feature in the critical be- 
havior which can be used to resolve the universality class. 
A striking feature of the RFIM is that the standard hy- 
perscaling relation between critical exponents is violated, 
and replaced by the modified relation 7 -I- 2/3 = ^{d — 9) 

urn 

Here, d = 3 is the spatial dimension; (3, 7, 



and are the critical exponents of the order parameter, 
susceptibility, and correlation length, respectively; while 
6* is a third (probably not independent ^3) exponent 
called the "violation of hyperscaling" exponent. For the 
pure (/3 « 0.326; 7 « 1.239; « 0.630 and diluted 
{/3 w 0.35; 7 « 1.34; v « 0.68 1]) Ising models hyperscal- 
ing is not violated, implying 9 — 0. In contrast, for the 
RFIM {/3 < 0.13; 7 « 1.7 - 1.9; v « 1.02 - 1.1 lH ^) 
one has 9^1. 

The aim of this Letter is two-fold. First, we demon- 
strate that hyperscaling is violated in quenched-annealed 
systems, using large-scale computer simulations and 
finite-size scaling. This fixes the critical behavior of 
quenched-annealed systems into the universality class of 
the RFIM, confirming the conjecture of de Gennes. Sec- 
ondly, we show that quenched-annealed systems also pave 
the way toward exciting new experiments. As we will dis- 
cuss, the considered quenched-annealed system is realized 
experimentally in a colloid-polymer mixture. Compared 
to aerogel, random-field effects should become easier to 
detect, due to the strong coupling between the fluid and 
the porous medium. 

The outline of this Letter is as follows. First, we intro- 
duce our quenched-annealed model. Next, simulations 
are used to show that hyperscaling is violated. In addi- 
tion, the critical exponents extracted from our data are 
shown to be compatible with those of the RFIM. Finally, 
we discuss how a quenched-annealed system can be real- 
ized experimentally in a colloid-polymer mixture. 

We consider colloid-polymer mixtures within the 
framework of the Asakura-Oosawa-Vrij (AOV) model 
[20I l2l| . The AOV model is known to reproduce ex- 
perimental observations remarkably well, including bulk 
phase-separation |22l |. interfacial properties 23], and 
gelation [iij. In this model, colloids and polymers are 
treated as effective spheres of diameter CTc and Up, re- 
spectively. The colloid-to-polymer size ratio is denoted 
as q = CTp/(Tc- Hard sphere interactions are assumed 
between colloid-colloid and colloid-polymer pairs, while 
polymer-polymer pairs can interpenetrate freely. The 
simulations are performed in the grand canonical ensem- 
ble, where the volume V and the respective (dimension- 
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less) fugacities, Zc and Zp, of colloids and polymers are 
fixed, while the number of particles inside V fluctuates 
(lengths are expressed in units of CTc)- Following con- 
vention, the polymer fugacity is expressed by the poly- 
mer reservoir packing fraction if^ = 'Kz^q'^ /Q. The colloid 
packing fraction is rjc — Tra^Nc/{6V), with A^c the num- 
ber of colloids in the system. In the absence of quenched 
disorder, the phase behavior of the AOV model is well un- 
derstood. For sufficiently large q, the AOV model phase 
separates into a colloid-rich phase (the colloidal liquid) 
and colloid-poor phase (the colloidal vapor), if rj^ exceeds 
a critical valu e n j [25j. The binodal exhibits an Ising 
critical point [2y|. The phase transition is driven by ij^, 
which thus plays a role similar to that of inverse temper- 
ature in liquid-vapor transitions of simple fluids. 

We now consider the AOV model inside a random 
porous medium. The medium is modeled as an equi- 
librium ideal-gas configuration of spheres of diameter 
(Tm = at fixed packing fraction rjM — 0.05 (in the ter- 
minology of Ref. I27L this resembles anticorrelated disor- 
der). We set q = 1.0, assuming hard-sphere interactions 
between colloid-matrix pairs, while the polymer-matrix 
interaction is left ideal. We use a cubic simulation box 
of edge L with periodic boundary conditions. We aim to 
measure the order parameter distribution [Ptlav Here, 
Pl = Pl{Vc) is the probability of observing a system 
with colloid packing fraction rye measured for one real- 
ization of the matrix, while [-Jav denotes an average over 
different matrix realizations. For each random matrix, 
a grand canonical Monte-Carlo simulation of the AOV 
model is performed, using a cluster move [2^. Colloids 
and polymers are inserted and removed from the simu- 
lation box at random, with the constraint that colloid- 
colloid, colloid-polymer, and colloid-matrix overlaps are 
forbidden. During the simulations, the (quenched) ma- 
trix particles remain fixed. The number of colloids in 
the simulation box fluctuates, and this is used to mea- 
sure [AF{Nc, Nc + l)]av, defined as the free energy dif- 
ference between the state with iVc and A^c + 1 colloids, 
averaged over typically 250 matrix realizations. By suc- 
cessively measuring the free energy difference, the to- 
tal averaged free energy [W]av as function of rjc results 
[2^ . The latter is related to the sought-for distribution 
[-PL]av oc e"'^^''^/'^'''^, with kB the Boltzmann constant 
and T the temperature. 

At two-phase coexistence, [Ptlav becomes double- 
peaked, where the peak at low (high) rjc reflects the col- 
loidal vapor (liquid) phase. The coexisting phase densi- 
ties follow from the average peak positions. Typical dis- 
tributions [Pijav are shown in the upper frame of Fig. ^ 
for several values of rj^. By recording the peak posi- 
tions as function of rj^, the binodal is found (see lower 
frame). Unusual behavior is revealed. For pure Ising 
critical behavior, well-separated peaks in the order pa- 
rameter distribution indicate that one is well away from 
the critical point, and inside the two-phase region of the 
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FIG. 1: top frame: Order parameter distribution [P^Jav 
for L — 12 and several values of rjp. lower frame: Binodal 
curves, obtained by reading-ofT the peak positions in [Pl]iiv, 
for several system-sizes L. The triangle, at 77c, cr ~ 0.070 and 
Vp.ci ~ 1.192, marks the location of the critical point in the 
thermodynamic limit. 
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FIG. 2: top frame: Cumulant analysis of the AOV model 
inside a porous medium at j^m = 0.05. Shown is Ui versus 
rjp for several system-sizes L. The vertical arrow marks our 
best estimate of 7?p,cr; the horizontal arrow defines Ui . lower 
frame: r;p,cr as function of 77M; the line serves to guide the 
eye. 



phase diagram. Finite size effects in the peak positions 
should then be small. In contrast, even though the peaks 
in [Pijav are well-separated for rjp « 1.17 and higher, the 
binodal continues to display a pronounced L-dependence. 
To locate the critical point, we have measured the L- 
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dependence of the cumulant Ui = [('Ti^)]av/[(|w|)]^^ with 
m = \ric — [(??c)]av|- At the critical point, the cumulant 
becomes system-size independent [23 ■ Plots of Ui ver- 
sus r]p, for several system-sizes L, are expected to show 
a common intersection point, leading to an estimate of 
77p The result is shown in the top frame of Fig. 13 
yielding rj^ ^j. — 1.192 ± 0.005. To determine the criti- 
cal colloid packing fraction 77c, cr, the quantity [(77c)]av(i) 
evaluated at rj^ ,,^, was linearly extrapolated in 1/L, yield- 
ing 77c, cr « 0.070. 

We thus find that the order parameter distribution 
at the critical point remains sharp, featuring two well- 
separated and non-overlapping peaks. Violation of hy- 
perscaling then follows from the L-dependence of the 
peak positions and the root-mean-square peak widths 
[sd Is ij . For the difference A between the liquid and va- 
por peak positions, finite size scaling predicts A oc L~^^'^. 
Similarly, for the width x of the vapor or liquid peak 
X oc L^f/'^^'^y^. The relative peak width thus becomes 
= x/A oc L'^ with uj = {'y /i/ — d)/2 + (3 /v. In case 
hyperscaling holds w = 0, and a finite relative width 
Wr > in the thermodynamic limit L — > 00 is retained. 
In contrast, when hyperscaling is violated and lo < w-c 
vanishes in the thermodynamic limit, leading to an order 
parameter distribution featuring two 5-peaks. Substitu- 
tion of the RFIM critical exponents indeed yields w < 0, 
implying — > for L — > 00, consistent with our ob- 
servations. Moreover, since the critical order parameter 
distribution tends to a sum of two (5-peaks, the value U* 
of the cumulant at the critical point (horizontal arrow in 
Fig. [SJ approaches the trivial value Ui — 1. Indeed, the 
cumulants of Fig. |21 intersect at a value close to one, con- 
sistent with RFIM critical behavior, but ruling out pure 
and diluted Ising critical behavior, where U* is signifi- 
cantly different from unity (see also Fig. O below). 

Having established that quenched-annealed systems 
show RFIM critical behavior, the experimental realiza- 
tion of such a model will be discussed. To this end, the 
crossover in critical behavior [s^, 113 , from pure Ising to 
RFIM, must be addressed. Defining t as the relative dis- 
tance from the critical point, the approach to the critical 
point, at < = 0, is characterized by two regimes: t < 
and t > tx, with tx the crossover temperature. For fluids 
in porous media, RFIM critical behavior is observed only 
when t < tx- For t > tx, the critical behavior is still dom- 
inated by the pure Ising model, and "effective" critical 
behavior is observed instead (in this case a combination 
of RFIM and pure Ising universality). Whether RFIM 
critical behavior can be observed depends crucially on 
tx ■ If tx is very small, precise temperature control is re- 
quired, which may be difhcult to realize experimentally. 
Note that tx is a non-universal quantity, dependent on 
the particle interactions, and, most importantly, on the 
packing fraction of the porous medium ?7m ■ By increas- 
ing TjM , random field effects are expected to become more 
pronounced, implying a larger tx- For silica aerogel, tx 
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FIG. 3: Evidence of the crossover in critical behavior, from 
pure Ising toward random-field Ising, as function of ?7m. The 
top frame shows U* as function of r]M ■ The lower frame shows 
the effective critical exponent ratio fi/v aa function of -qm- The 
dashed lines serve to guide the eye; horizontal lines show pure 
Ising and RFIM values. 



is clearly very small, since the measured critical expo- 
nents do not differ much from those of the pure Ising 
model In contrast, our results for the quenched- 

annealed model show pronounced RFIM behavior, indi- 
cating a much larger tx- For experimental applications, 
it then becomes relevant to know which values of 77M are 
required, in order to enable measurements in the regime 
t < tx and observe RFIM behavior. 

Shown in the lower frame of Fig. [3 is r7p cr as func- 
tion of r7M 01 • Defining the analogue of Eq.(^ as 

S = \vl„{M) - vlAP)\/vi,cr{P)^ ^itl^ vlcAM) the 
critical value of 77p in the presence of the random ma- 
trix, and ?7pcr(^) fti^ corresponding value in the pure 
(77M = 0) system, S is found to increase from S « 0.03 
(77M = 0.005) up to <5 « 0.5 (77M = 0.05). This confirms 
our expectation that, by increasing 77M, random-field ef- 
fects grow stronger, and tx becomes larger. Evidence for 
the crossover in critical behavior is presented in the top 
frame of Fig. O Shown is U* as function of 77m. Recall 
that Ui is defined as the value of the cumulant at the 
critical point, see the top frame of FigJEJ The horizontal 



lines correspond to Ui 



1.2391 |2J| of the pure Ising 



model, and the (exact) RFIM value 11^ = 1. The figure 
strikingly illustrates effective critical behavior, between 
that of the pure Ising model and the RFIM, with a pro- 
nounced drift toward the latter, as t^m increases. Addi- 
tional confirmation of the crossover is obtained from the 
critical exponent ratio P/v. Shown in the lower frame 
of Fig. |3is l3/iy as function of 77M. The upper horizontal 
line shows the pure Ising value; the lower line is an upper 
bound for the corresponding RFIM value [l8lll9| . Again, 
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a clear drift toward the RFIM value is observed. 

Fig- El then provides a clear indication which value of 
77m to use in an experiment. The relative distance from 
the critical point that can be reached in simulations is 
nowadays t « 10""^ (sll; similar precision is also achieved 
in experimental colloid-polymer systems @|. Although 
this precision is rather low compared to what is achieved 
in atomic fluids, Fig.Onevertheless shows pronounced de- 
viations from pure Ising behavior already at 77M = 0.015, 
with the crossover to RFIM being nearly completed at 
r]M = 0.05. This suggests 77M > 0.05 as the optimal 
regime for experiments. Such packing fractions are sur- 
prisingly low, and easily realized in aggregated colloidal 
rods or spheres . There may even be the exciting 
possibility to generate the porous medium by optically 
trapping some of the colloidal particles [s^. An addi- 
tional advantage is that colloidal particles, due to their 
mesoscopic size, allow for very detailed investigations of 
fluid phase behavior. By using confocal microscopy, in- 
dividual particles can be visualized and tracked directly 
in real-space ^ . This has already enabled particle-level 
investigations of interface fluctuations [ij and bulk crit- 
ical behavior [s^. The present results indicate that the 
experimental verification of random-field behavior is fea- 
sible in colloid-polymer mixtures. 

In summary, we have used large-scale Monte- Carlo 
simulations to resolve the universality class of the 
quenched-annealed model. The universality class was 
shown to be that of the random-field Ising model, as was 
evident from the violation of hyperscaling, and the be- 
havior of (effective) critical quantities. This confirms the 
conjecture of de Gennes. In addition, we have demon- 
strated the potential of colloid-polymer mixtures in the 
experimental detection of random-field critical behavior, 
providing a valuable alternative over aerogel-based sys- 
tems. 
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